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SUMMARY 



A new method is presented -for obtaining tho velocit;? 1 
potential of tho flow about a solid of revolution moving 
uniformly in the direction of its axis of symmetry in a 
fluid otherwise at rest* This method is "based essentially 
oh the fact that the form of the differential equation for 
the velocity potential is invariant with regard to conformal 
transformations in a meridian plane, 3y means of the con- 
formal transformation -of the meridian profile into a circle 
a system of orthogonal curvilinear coordinates is obtained > 
the main feature of v/hich is that one of the coordinate 
lines is the meridian profile itself. The use of this type 
of coordinate system yields a simple expression of the 
boundary condition at the surface of the solid and leads to 
a rational process of iteration for the solution of the 
differential equation for the velocity potential, . It is 
shown that the velocity potential for an arbitrary body 
of revolution may be expressed in terms of universal func- 
tions which, although not normal, are obtainable by means 
of simple quadratures. 

The general results are applied to a body of revolu- 
tion obtained by revolving a symmetrical J.oukowski profile 
about its axis of symmetry, A numerical example further 
serves to Illustrate the theory. 



INTRODUCTION 



The simplest case of a three-dimensional fluid motion 
occurs when a body of revolution moves with a constant 
velocity in the direction. of its axis of symmetry. In 
this case,. the motion is the same in any plane passing 
through the axis of symmetry and, in this respect, presents 



2 



some analogy with a twp-dimens ton'al motion. Thus, a stream 
function is defined "by Beans of the equation of continuity; 
and the condition for irrotational motion yields a velocity 
potential. The stream ' function', and the velocity potential, 
however, are not int erchangeable in three-dimensional flows 
in the same way as '-are t he corresponding quant it i es in two~ 
dimensional irrotational motions. The reason for this 
dif f erence x i s "that, although the differential equation for 
'•the velocity potential is Laplace *-s equation, the equation 
for the stream function is not Laplace's and therefore 
the two functions- cannot ; "be combined' tc give an analytic 
function of a single complex variable. It follows that 
the elegant and powerful methods of the complex variable 
are not obviously • applicable 'and' the- calculation of the 
pot ent ial flow, past a body of revolution has-, of necessity, 
developed alorig other lines. 1 " 

The method of calculating the flow past a body of 
revolution most often referred to- was suggested by Hankine 
and developed by von Zdvmdh (reference 1) and others. 
For axial flow, the axis of the body is coveredby.a 
continuous distribution of sources' and sinks in such a 
way that the closed stream surface found by the super-' 
position of the flow induced by the sources and sinks on 
the parallel flow coincides with the surface "of the solid. 
For transverse 'flow', the. axis of the body is covered by a 
continuous distribution of doublets.- The superposition of 
the flow due to. the doublets on- the parallel flow. then- 
yields the surface of the solid as a stream surface. Both 
problems lead to integral equations that von IC^rmdn solved 
by an approximate method. As von.K^rm^n pointed-out, how- 
ever, the exact replacement of the body by a distribution 
of singularities along the axis of symmetry is possible 
only when the analytical cont inuat ion of the- potential 
function, . free from, singularities in the space outside the 
body, can be extended to -the axis of symmetry without en- 
countering singular .points . (See reference 1, p f 27.) 
Inasmuch, .as the question of when this analytical contin- 
uation is or is not possible has never been answered, 
numorical calculat ions ' may .lead to incorrect results, 
particularly in the case of a body with a rather blunt 
nose. , 

•In reference 2 an attempt was made to calculate the 
potential flow past a body' of revolution according to the 
methods of potential theory. In- that -paper , Laplace *s 
equation for the' velocity ■ potent ial is expressed in terms' 
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of elllpt ; ic-;C*3rlindrical coordinates and is. solved in 
conjunction, wi.th the. appropriate "boundary conditions for 
• axial, and .transverse f lows. In the case, of an ellipsoid 
or a hyperboioid- of revolution -the solution is obtained 
co in a closed, form, .."both shapes "being member s:- of the family 

in of orthogonal coordinate surfaces belonging to a system 

►i of elllpt ic-cylindr ical coordinates * Por : any other, solid 

o.f revolution, howev.er, the method leads : - ..to t.wo sots of 
..linear equatians , each set. having an infinit-o number of. 
equations and an infinite number of unknown coefficients 
f or .the determination of the veloc i.ty .po t-ent i-ais. for the 
axial and, t.-ransvor $e flows. .... j. -.• 

In the present paper, a new method is presented for 
calculating the' potent i-al f low -past -.an. arbi trary body of 
revolution. Only the case of axial flow is discussed but 
the methQ:d is equally, applicable to the- case of transverse 
flow. She method .is based on the discover:/ that, by the 
proper choice, of .the system of coordinates to be. used for 
a given : bod3" of revolution,- the: solution of the potential- 
flow problem' can be o-otain^ed by mean-.s of- quadrat.ur.es. - 
Coordinate systc-ms of th is.- nature already, exist in the 
literature., for example if the: solid is an; ellipsoid 
or a. hyperboioid of revolution, the. coordinate system 
used is an ;c 11 ipt ic^cy lin-dr.ical one*- - These . coordinate 
systems: possess in common the .property that one of - their 
co-ordinate surfaces . i:s the boundary of the - solid itself.. 
It will be shown in this paper that, by means of the con- 
formal transformation of the. meridian profile into a circle, 
a system of orthogonal curvilinear coordinates can be de- 
fined such that one of the coordinate lines is the meridian 
prof lie. of.;-. the solid of re volut ion*. . Fur the r more , it will 
be shown that, by the use of ' this typo of coordinate system, 
the potent inl-»f low problem for an arbitral body of ..rev.olu*- 
tion can bo .solved by means of elementary integrations. 

In comparison of the method of this paper with the 
..method of reference 2, it. is to be recalled that, in refer- 
ence 2, Laplace ' s equation for the Velocity potential was 
expressed in elliptic-cylindrical coordinates. The general 
solution involves normal functions:, namely, Legcndrc func-* 
t toffs of the' first -and ■ second kinds-. The -'boundary condition, 
however, for an arbitrary body of revolution' - involves the 
two independent variables of the problem and yields, as 
mentioned before, an infinite set of linear equations with 
v an infinite number of unknown coeffic icnt s. On the other 
hand, the method of this .paper..: ut ilizes a ■ different " co- 
ordinate system for each body of revolution. The main 
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feature of .these conf ormal orthogonal curvilinear co- 
ordinates is that one of the coordinate lines coincides 
with- the meridian profile of the "body of revolution. This 
fact leads to. an expression of the "boundary condition at 
the surface involving only one independent variable . 
Although the general solution of the differential equation 
•for the velocity poten't'ial . has not been obtained in -terms 
of normal functions, a me,thod of iteration has been de- 
vised that involves only elementary integrations. This* solution, 
satisfying the boundary conditions, can- then be expressed 
in- terms of universal functions.' Although these functions 
are not normal, they need be determined but once. 



MATHEMATICAL "DEVE-LOPMEHTS 

-Equation for the Velocity Potential 

The axis of symmetry of the body is denoted by x 
and the. position of a point in a' meridian plane is fixed 
by .the Cart osian ' co-ordinates (x, Q>) (fig. 1). Then if 
c[ x and qjrj are, respectively, the components of the 
fluid velocity in the directions of the - x. and U5 
axes-, the equation of continuity is obtained by equating 
to zero the flow out. of the. annular space obtained by 
revolving a, small rectangle dxduJ around the axis of 
symmetry. ' Thus ,. the total flow outward in the direction 

of the x axis. is. J^(2ttu5 q. : . : duJ ) dx and in. the direction 

6x . 

of the . w axis is ~r ( 2ttu5 Omdx) dtp. The equation of 

ocu ^ 

continuity is therefore 



~ + ^ („- ftjB )-= 0 . (1) 



Since the flow is symmetrical about 'the x axis, • 
the vorticity is 



Ail 

cix 



oUJ 



and if, fxirther, the motion is irrotational , then 
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A velocity potential 



can therefore "be so defined th 



and the equation of continuity (l) becomes 

4i. f a) ] + . { -£ = 0 
ox V ox; da) \ ouj/ 

Consider now the conjugate complex variables 
z = x + iu) and z* = x iu>. Then, symbolically, 



ox 



1 .Su) 



?Fz §x + ' x 6o» 



and it can "be easily verified that 

Thus, the vanishing of the real part of -j- (cJ ~ = J is 
equivalent to equation (3 ? }>. , v ® Z/ 

Consider further the conformal transformation 



z = f ( £) ' vhere J = | + in 
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Then 



» (, m ♦ » (, ||Y . S .P. 4 it » ( 5 SI |£) (4) 

Since z, z, or t » ^ • may he looked tipon as. independent 
variables, it follows that the right— hand side of equation 
(4) is 



Mow, the -product is a real quantity. The van— 

dz dz 

ishing of the i\eal part of ~~ (jS ^f) is therefore 
equivalent to equation (3). Hence, 



d| V- hiJ hi] \ bn 

where uJ is now a function of |, n obtained from the 
conf ormal tr ansf ormation z = £(£)•* 

*It has "been pointed; out "by Dr. Theodore Theodorsen 
that equation (5) can.be obtained directly, as follows: 

The vector form of equation (3) in a mer idian- plane . 
is s imply 

div (u) gr ad 0 ) = 0 

If a set of orthogonal curvilinear coordinates \ and rj 
is introduced, the expression for div (cu graft 0) becomes 

1 V Y m h s W\ ; a / h x &0n" 
si \* 17 ay; + h~ 

where the elements o.f length along the variable and r\ - 
variable coordinate line's are, ro spect ively , h x d^ and h 2 dr). 

If, however, the transformation from the rectangular 
Cartesian coordinates to the orthogonal curvilinear 

coordinates (|,r)) is also c onf ormal then h x »h 2 ; because 
for a conf ormal transformation the magnification at 'any 
point in the plane . is independent of the 'direction through 
the point. Thus, with h x *s.h at the foregoing expression 
placed equal to zero yields equation (5) # 
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Equation (5) rather ' tftan ' equat ion (3) will be con- 
sidered to "be. the fundamental ..differ ent ial "equation 
governing axisymmetr ical mot i p&s ' of / a .per f eet iricompr es s— 
'ible fluid> It is more general th^'n equation (3) . inV the 
sense that' the independent ..variables (|, tj") denote any set 
of or thogonal curvil inear coordinates' obtained "by "means of 
a conf oraal transformation z = f(£')» In addit ion,, • it s 
form is invariant with regard to conformal -transformations 
and, in this sense, does not complicate the original 
equation (3). In the following section it will be shown 
that, by means of the conformal transformation of the 
meridian profile ox an arbitrary body of revolution into 
a circle, an orthogonal curvilinear system of coordinates 
(|, n) can be so defined that the coordinate line iy = 0 
is the meridian profile;, itself ... -In .subsequent sections 
it will be seen that the- use of this type of coordinate 
system leads to a simple' expression: of the boundary con- 
dition and, consequently, to an iteration process for the 
solution of equation (5) for an arbitrary .body .of ; rpv ; olut ion. 

Conformal Trans form at ion .and Orthogonal Coordinates 

It is well known, that a unique •conforma! tra^sf ord- 
ination exists which maps the region exter-nal to a given 
boundary in the z plane into the :• region . cxt ernal to a 
circle in the. Z piano wit h,. it s ' cent or at the origin, 

' feuch -that the r egions-at- infinity of the two-. planes 
correspond. Ehe. f unc t ion r opr ps ent ing" thi s conformal 

'transformation can be' devclbpbd' in ' the 'region external to 
the..? circle in a. convergent, series , of thp type, . , 

a i a 3 a 3 . ' . v * ' ■ 

z = Z + c x + -2- js * * \ (6) 

where the coefficients c x , a. x , a 2f - a 3 , . • • are, in. 
general, complex quantities. -.In this paper, the 
meridian profile is symmetrical with respect to the 
axis of revolution and these coefficients, therefore, 
are all real . The constants a a , a s , a 3 , . • .depend 
only on the shape of the nveridian profile in the z 
plane. 

' The method'-, of Theodorsen and Garricfe, described in 
• reference 3, is par t icularly • well, adapted for the pur- 
pose of determining transformation (6'). "It ' is ; shown in 
reference 3 that, with the proper choice' of axes and 
origin, the Joukowski transformation 
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maps a closed "boundary in the z plane into a nearly 
circular .curve in the Z 1 plane, The mapping- of this 
nearly circular "boundary , v/ith its center at the origin ' 
of the Z plane, is then completed "by sieans of the 
transformation 

Zi = Ze f(Z) ( 8 ) 



where 

«.> - j '51 

n=l 

On elimination of Z 1 between equations (7) and ( 8 )\ 
equation (6) follows ^ where 

1 

H * c 2 + ■§ + a 2 

1 „ • 

a 2 = c 3 + c 2 C! + 5 Cl 3 - Cl aa 



The valties for c n and the-refore for a n • can actually 
"bo determined as the Fourier coefficients of a certain 
\|/(<p) curve (refer once 3). 

If the radius of the circle in the Z plane is 
denoted "by E r the* coordinates • of ^ any point '.on this 
circle can "be expressed as 



X = R cos i 
Y = -H sin I 



where, as £ increases from 0 to 2rr , the point describes 
the circle in a clockwise sense, thus leaving. the external 
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region to the left of the direction in which the boundary 
is traversed* Any point of the circle can also "be 
expressed as 

Z = X + iY = Re"~* ^ 
If t - t + IT] , the transformation. 

Z = Re'^ 

yields the circle of radius R for ri = 0 and, further- 
more,, n = o° corresponds to the region at infinity of 
the Z plane. Upon substitution of this expression ' ■ 
for Z in equation (6), it follc\/s that 

it R = T£3 



When this equation is separated into its real and 
imaginary parts, it is seen that • 

x = c 1 + Re* 1 cos | + —e""^ cos g + S| e ~ Sr l cos 2| ^ 



+ 3 e -3t1 cos 3|+ . . 



a) = -Re 11 sin g + -g-e""^ sin | + g^e-STj sin 2g 



where the "quantities c lf a l , a 2 , a 3 , . . . are deter- 
mined according to the method of .reference 3. 

Rather than as a conformal transformation of a 
plane z into a plane Z., equations (10) are to be 
looked upon as the equations of transformation from the 
rectangular Cartesian coordinates (x, uJ) into the orthog- 
onal curvilinear coordinates ( g ,n K Furthermore, the 
coordinate line r) = 0 is the .profile itself; that is, 
when r\ = 0 , equations (10) yield the parametric 



> do) 



J 
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equations of tho profile in the piano z. It will bo 
scon in tho 'following.*. discussion that the use of this 
type of coordinate system leads to a sirnplb form for 
the boundary condition at tho solid surface 



Boundary Conditions 

If a body, of revolution moves with a velocity U 
in the direction of its axis of symmetry, the normal 
velocities of tho body and the fluid in : contact with it 
are the same. Thus, at the surf ace of,. .tho mov tng body , 
the boundary condition expressed in the .conf orma l ortho- 
gonal coordinates (|, n:) is simply 



Tl=0 



ri=0 



since tho coordinate linos along wh'ich ti varios and 
I remains constant arc normal to tho boundary ti =' 0. 
According .to tho first of equations (10) this boundary 
condition can be written as , 



-(*$) =U (Woe I- Si cosl- 



2a s 



cos 2| ---^cos 3|- # .Mll) 
H3 J 



Furthermore, as the fluid at inf inity , originally . at rest, 
remains undisturbed by tho' motion of tho body, the bound- 
ary condition there is, simply, 



(8).. 



= 0 



and 



> 



.= 0 



(12) 



It is. remarked that the simple form of the. boundary 
condition (ll) has' been attained on account of the 1 intro- 
duction of conf ormal or thagon'al coordinates (£ f tO'so 



11 



defined that one of the Coordinate lines ''ti' = 0 is 
the meridian profile itself, I-t"-will # be seen in the 
following discussion that this choice of coordinate 
system lead> to a. process of .iteration for , the solution 
of the ftindamehtal differ ent ial equation(5) involving 
only simple quadratures. 

Solution of the Fundamental Differential 
Equation by Iteration - 

If the right— hand side of equat ion: ( 10 ).. f or to ' 
is substituted into equation (5), it follows that 

n-2 u ■> 




(13) 



With regard to an iteration method it is desirable that 
the initial step in the process bB obtained in a closed 
form. Thus, equation (13) has been so arranged that the 
solution of the left— hand side placed equal to zero can 
be obtained in a closed foriii. This- init ial solution can 
then be utilized as the starting point of an iteration 
process. Before a. detailed/ description of the iteration 
method is given, hoxvever,' it is first necessary to in— . 
troduce several new parameters. Thus, the coefficient 
a x is replaced by a 3 e 2 ° and the radius E by ae a 
where a,, as in equat-ion ( 7 ) , serves merely to preserve 
dimensions. Then oquation ( 13 ) can be written as follows i 
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sinh ( n+a-P ) sin | (*-J^ + T~4Y + sinh (n+a-0) cos % ^$ 

\df; s otj2/ of; 

+ cosh U+a-P) = Y h n e "( n+1 ) (n +a -P) 

n=i 

X £sin (n+l) 1 + ~|) + (n+l).'cos (n+l) t || 



- (n+l) sin (n+l) % ~ | .(14) 



&0j 
J 



'• &n+l . 

where d =• — ■ a -— , a nondimensi onal quantity, and 

2(ae p ) n+<! . 

the "boundary condition ( 11 ) . "bocomos 

- (M) 2aUe P sinh (a-P) cos' | 

- 2D 1 e^ 2 ( a "P) cos 31) 2 e~3(ar-p) C os 3|- # ...J ( 15) 



The method of solving equation (14) is based upon the 
following cons id'er at ions : 

It is assumed that the velocity potential 0 'can 
he developed in a series of terms each, of which is homo- 
geneous with respect to .the - indices of the coefficients 
b n ; that is , 



CO 



^0 

where , for example, ' 

0 O involves none of the coefficients h n 

0 X contains the coefficient h^ as a factor 

0 2 is the sum .of two terms having, respectively, h x 2 
and t> 2 as factors 
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0 is tho sum of three terms having, respectively, 
3 t,^, "b 1 'b 2f and b 3 as factors and so on for 

the higher-order terms 

Equation (16) for 0 is now substituted into equation 
(14) and the terms involving the coefficients b n and 
their products to the same dogrce in the indices are 
equated. The same process is applied to the boundary 
condition (15), In this manner, equation (14) is re- 
placed by a sot of part ial dif f or ent ial equations which, 
in conjunction with the boundary conditions, can be 
solved rigorously for 0 O , #x> 03. • • > VhGn tko^ oper- 
ations just described are performed, the differential 
equations with the accompanying boundary conditions for 
tho first threo functions. 0 O • 01 » and 02 arc as follows 



1 
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< Determination of fo 

In order to solve equation (17) for 0. o ,. it is. 
convenient to introduce a new set of independent variables 
p and X , wh er e 

p, = cos i and X = cosh (Ti+ctr-p) 
By use of the symbolic expressions _^ 

Tr = -sin % -r- 



and 



*1 = -oos I + sin*! ^ 
d| dp, op, 



= sinh (n+a-p) 
OT) ^\ 



> 



= cosh (ti+co-P) — + sinh (n+a-p) -r-g 
on °* 

it follows that equation (17) can be written as 



(20) 



dp, 



dX 



= 0 



(21) 



If a solution of the form 0 O = '?(p.) &(X) is assumed, 
equation •( 21 ) becomes 



Since' the left-hand side of this equation contains p, 
but not X and the right-hand side contains X but not 
u, F(u) and G(X) must be such that each side is a con- 
stant. If, furthermore, the constant is chosen to be 
-n(n+l), where n is a positive integer, then 
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(22) 



A. ["(l-\ 3 ) + n(n+l)0 = 0 



cLXi ! ' d\ 



Each of these equations is a Legendre differential 
equation and, "being of the second order, possesses two 
independent particular solutions: namely, the Legendre 
functions of the first and second kinds, P^ and Q n> 
respectively,' According to the "boundary condition at 
the surface, the form to be taken for is' Pida). 

Correspondingly, according to the "boundary condition 
at infinity, the form to he taken for &(\) is ^Q^Ot 
since the Qn(X) and their derivatives' vanish for \ = <»« 
It follov/s that the solution of equation (2l) is 



The arbitrary constant A x is determined by .the boundary 
condition at the surface-; namely, 



(23) 




Thus , 



2aUe p 



wh er e 




and 



\ Q =• cosh (a— 3) 
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It will he necessary for the determination of 0 X 
and fl5 3 to have the first few orders of Pn(M-) ancl - 
Q n (X) stated explicitly* The following are expressions 
for the first four harmonics:' 



P 0 U> = 1 



Pi(lJ-). = P- = cos i 



PsCjx) = |C3n S - 1) = .| cos 2 | + | 
.'^(P-) = i( 5 U< 3 - 3p>) 3 f cos 3't + § cos i 



> 



q 0 (x) .- | log ^ 

r 



' Q 2 U) = |(3X 2 r. 1) log »™ - | X , 
. Q 3 U) .-. J(5\ 3 - 3X) log - | X 3 + | 

Determination of 0 X 

When the independent variables' |j, * and X are 
introduced and equations (21) and (23) are used, it 
follows that equation (18) for 0 X can "be v/ritton as 



J 



(34) 



3 /■ y-\» - i ^ v . J 

- 2|"<il(X) + (X + l) Qx'CXjJ P 2 (|x)j (25) 
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The right-hand side of this equation requires that 



0x. --.F 0 (\.) ? 0 U> * ?. 3 (X.) F 2 (>) 



(26) 



If :ihis expression., for . ; 0 X •; is sutist itut ed . into equation 
(25) and th e ':. coef f ic lent s : of .'P 0 (fx) . and . P 2 ( jj,) .are equated 
oh' both sides, of the ..r e suit ing equation,- it follows that • 



d 

d\ 



dF 0 (\) 



2 
3 



b^ 



dX, 



f 7>? -i . v 



(x+2.y ^ 1 ( ^ )-2Qi.(x ) 



(27) 



and 

_d_ 
dX 



,-d\ J 



3 ■ j\z -l . l> / . ■ - J 



(28) 



With the right— hand sides t ak en " equal ' t o zero',"' : th'ese two 
equations reduce, respectively, t o* Legendr e 1 s equation for 
the zero- and .second-order functions. Then, according to 
the theory of linear d iff er. en t ial. equatlp.ns ..of the. second 
'order, the s ol\it ions of cquat ions .( 27 ) .and ( 2,8' ) aror , # 



1 Q M * B 0 Q 0 (\) - 



d\ 



/R c (\)cl\ 



(29) 



and ■ 
2? 2 (X) = B 2 Q 2 (X)-(3X s -l) 



d\ 



(X S -1)(3XS_!)' 



/(3X 2 <-l)R 2 (X)dX (30) 



f 
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where R 0 (\) and H 2 (X) are, respectively, the right— 
hand sides of equations ("27) and ( 28 ) . It. is noted that 
each of those solutions contains only ono arbitrary con- 
stant, 3 0 and B 2f bceauso the other two independent 
solutions P 0 (X) and' P 2 (X), respectively, of the homo- 
geneous forms of equations (27) and (28) do not satisfy 
the boundary condition at infinity and therefore the 
arbitrary constants associated with them are taken equal 
to zero, . 

The boundary conditions for 0 X lead to the following 
boundary conditions for F 0 (?0 and P 2 (X.)ef 



(*o\ = V — 'J * 0 



(31) 



and 



16 

-g- aUe 



2'\=co \cUA_o, 



p( x o-y^o 3 -i) 3 ; 



(32) 



The int egrat ions required in equations (29) and 
(30) are straightf orward and need not be performed here # 
Furthermore, by use of the well— known recurrence formulas 
for the Legendre functions, it can be shown that 



and 



- I D X A X || [2Q 0 (X)+Q 3 (X) J + yX B -l Q 0 ! (^) * 2 ^l(^)]j , (34) 
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where, "by means of the boundary conditions, equations 
(31) and (52), 

1 



B Q = 0 



and 



B 2 = GbiAj 



Q 3 '(x 0 ) 



(35) 



It may "be easily verified that these expressions for. 
F' 0 (\). and F 2 (X) satisfy the condition that they and 
their derivatives vanish for \=oo. 



Determination of 0 3 

When. the variables u, 'and \ are introduced into 
equations (21) and (35), it follows that equation (19) 
for 0 S can be written as 



A - 1 



Mr( 1-}1 2 ) 







6X J 










- M 






dp, 




60o 





Of* 



^ + - 2 yPii)Mo](36) 



If the expressions for 0 ? and p x given, respectively, 
by equations (23) and (36) are substituted in equation 
(36), it follows that: 
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The r ight— hand side of equation (37) requires that 

When this expression for 0 2 is subst itut ed .int o equat ion 
(37) and the coefficients of Pi((x) and P 3 ( X ) are equated 
on both sides of the resulting equation, the following 
equations for ^(X) and ]? 3 ( X.) are obtained: 

^ [(l-\ S ^^^] + S^U) ^(X), (39) 



and 



_d_ 

dX 



d M X >' + 12F 3 (X) = H S (X) (40) 



It is. evident that the homogeneous forms of these 
equations are, respect ively, Legendre 1 s equations for 
the first— and third— order functions. According, then, 
to th.c theory of linear differential equations of the 
second order, the soiut ion«s: of equations (39) and (40) 
ar c '*■' ■'■ 

*i(X)-= B^iCX) - /H x (X)dX (41) 



and 

g 3 ( X ) = B 3 Q 3 ( X ) - ( 5X3 _ 3X _ . ( d ^ 3 _ 3x } 3 /( 5X» -3X )E 3 ( X ) dX 

'# ( 42 ) 

The arbitrary constants B x and B 3 are determined by' 
the boundary condition at the surface of the moving body. 
Thus, the boundary condition for 0 2 leads to the follow- 
ing conditions for ^i(X) and I3 ( X ) : 



0 
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and 




Then,, if the integration^ : re:quiire:d in equations. (41)^' 
and (42) are made- and; Vtefe 'recurrence; f o rmul as.? -f o r : . ' t he! 
i^g-endre functions' are Used,': it can be •shown .{hat 
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' f?^!^ 4 ^! <M +§^5 (X ) -VX 2 = 1 JOg (X ) *Qo (X ) - gQi'CX)]^ . (45) , 



and 



P 3 (X ) • = B 3 Q 3 (X 5 - 0% ^ } ♦§% > -h^l [d^ (X } + 4Qg (X jjj 

- ^^{^l^)- ^aJ^X^Tl j§Q 2 (X)- |Q 0 (X)+%M'Mjj (46) 



- for^JsQgfXoH^j!' (s^f^ ?f ^o>-|v^^pq 1 {X oK Q^ (X 0 )=f (xj 



and 



^^^(X^^!^^ (49) 
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It may "be easily verified that the expressions for 
? X (X) and (M and their derivatives vanish for 
A = oo • 

It is clear that the iteration can he continued 
to obtain further members of the sequence of )6 n . The 
order n to which the process must be carried depends 
on the magnitude of the coefficients b n of the conformal 
transformation, For example, it may be that one of the 
higher coefficients, say b m , is significant. It is 
obvious, then, that the iteration must be continued far 
enough to include^ at least the corresponding term 
Although the 0 n f s.do not form a set -of normal functions 
(obtained independently of one another) t they are never- 
theless uniquely determined by means of simple quadratures; 
that is, they constitute a set of universal functions in 
the sense that the labor necessary for their determination 
is performed bu-. once and need not. be repeated for, each 
particular body of revolution. 



Velocity and Pressure 
According to Bernoulli's theorem, 




(49) 



where > 
p pressure anywhere in fluid . 

Poo pressure 5,n region at infinity undisturbed by 
motion of solid 

p density of fluid r. 

<i r velocity of fluid relative to boundary of solid 



It is recalled that in this paper the expression for the 
velocity potential, 0 has been derived by considering 
the coordinate axes as rigidly attached to the body, and 
therefore moving with a velocity U " in the positive 
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direction of the x axis. The velocity vector 

q = - grid 0 represents the velocity of the fluid 
relative to the undisturbed fluid at infinity. The 
vector Telocity q^^ relative- to the moving axes, is . 
therefore given "by 



or \ / 

t r = i(q x " U) + J q£ 

where ' i and j are, respectively, unit vectors along 
the positive directions of the x and <5 axes. The 
magnitude of the relative velocity is therefore obtained 
from 



q r 2 = q/ + 2 Uq x + U2 (50) 

where 

<30 , S0 

It is. recalled, however, that the velocity potential 0 
is derived as a function of the independent variables £ 
and T) • In order then to determine q r it is necessary 
to express the velocity components q x and in terms 

of the quantities and These expressions can be 

or ; -an 

derived in the following way: 
The relations 

z * x + i<5 and z « x - 125 



,can be considered a3 equations of t reins formation from the 
coordinates (x,o) to the independent coordinates (z 9 z) 9 
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It follows then-, symbolically that 



2 JL = JL - i JL 

5z ■ dx . 



ana 



Ther ef ore 



2 4- . 4- + i.-4. 

d"z ox das 



a + i q_/= - (± + i -L^ 0 = J.2 ^ 



whore 0 is now a function of the indopcndon t variables 
z and "z. 

Iho conformal transformation . z = f (.£ ) then gives 

• ■**■+ /.to. = - 2 . §f« ...... 



or 



„ <L d.z 

~ ' J g d£ d£ *' 



(51) 



where 



Symbolically , 



2 d % d z 

J ".IT it 



2 A = A 



i iL 

Sri 



and therefore 



= ( TT + — J + i h 

\ ot on / 



dl \ of OT| 



) 
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Then, by means of the Cauchy— Hiemann relations 



and 



it follows that 



or 



d z 6x 



_6x 
8ti. 



oar 



A 2 _ bm . 
d£ " b n 1 6 | 

Equation (5l) therefore becomes 



+ iq- 



U) 



" St) A; 



i r=Atf + .1 . 

Ot OT| 



(53) 



where 



" ■ (if r ♦ (n y - (f r* (f t ■ 

When the real and imaginary parts on both sides of 
equation (52) are equated, it follows that 



* x ; Vat al. a n 



and 



> 



1 .da? •■ &ff duaN 

" r " ~ j 2 \ai "&e ail an / 



(53) 



By means of these, equations , -q r ; 3 " and therefore 'the pressure 
•coefficient 



— - — ^ given, respectively, by equations (50) 
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and (49) can. be. obtained .as functions of the independent 
variables \ and n 9 ' 

APPLICATION OF THEORETICAL RESULTS 
Symme trical J oukowski Shape 

General e xpressions The theoretical' results of this 
paper will now be applied to the case of a body of revolu- 
tion whose meridian section' ts a symmetrical Joukowski 
profile with both rounded nose and tail. This example 
is sufficiently complicated to illustrate the principles 
and usefulness of the method. 

It is .well known : .that by means of tne mapping function 

z = z * + £7 (54) 

the circle of radius a, with its center at the origin 
of tho H) plane, is transformed into the line segment 
extending from 2 = 2a to z - 2a in th.e z plane; 
and the circle of radius 1, a(l +' c 1 + €2) with it s 
center at Z t = e x a is transformed int'o a symmetrical 
Joukowski profile with rounded nose and tail in the z " 
plane, (figt 2), When Z x is replaced by Z + c 1 a t 
equat ion ( 54) becomes 

z = Z + V x a + a3 '• 



or 



S 3 . 4 5 

a 1 a' ■ . 2 a 3 a 



z = Z + € X a +■ y-'ef^ + e-i 2 "gs ^ e i 3 . g4 + • • • ( 55 ) 

A comparison of this equation with equation (6) shows that 
c x = ^a, a x = a 2 , a 2 = rCia^, a 3 = e^a 4 , a 4 - — e^a 5 ^ 

A 

Also',' the radius of the circle in the Z plane being 
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denoted by ae a t it follows that 

e a = i + €l + €2 (55) 

According to definition 

2 b H ' a m+t • ■'* ■' " 

© p = — and b ra •= — a ' g 

■ a 2 m 2(aeP) m t^ ; ' . ,-. 

wher e 

a > 1 

As a lr ~ a 2 .and . a m+1 ~ ( T l A 1 m a D1+? ,. it follows there- 
fore that ' . ; ! .' : 

p = 0 and b a = -(-1) 4 (57) 
The expression for the velocity potential then becomes 

0 = 0 O + 0i + 0a+ • • • (58) 
whero, from equation (23) 

0o = A^UJC^U) (59) 

with; 

2aU 

= - 

and 

\ 0 = cosh a - j (l+€i+€ s + 1 +t 1+ £ 2 ) 
Prom equation (26) 

01 = F 0 (X) P-(n) + P 2 (X).'P s (|x) (60) 
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where 

» 0 (X) = ^i A i{| [2Q 0 (X). + Q 2 (X)| . 

+ fi7~^i [vOO - Qi(x)]} 

I 8 (X) =B 8 Q a (X) +|c 1 A x |i [3Q 0 (X) + Q s (X)j 
+ V'x 3 - 1 [V < x ) + 2^1 < x )| 

and 

BQx'^i) + syt/'tTT ('q. 0 .(x 0 )V 8 



b. 



o > Yx. ) 
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Velocity and pressure From equation together with 

the transformations 



-i t 

Z' = «j_a + Z and Z = ae a e~ 



it follows that 



= a(€ 1+e 71+d cos i) fl + 



M62) 



and 



© = -ae^ +a 3in£ 



Then 



e 2 2(n*o) 
*1 e 



and 



T)+<* - , 2(ti+ct)j 2 , 

1 cos§ ■+ e 



r 2 + e 2(n + c)j e ^a co3| + 2€ie 2(^a) 



€, 2 + 2£,e^ a eosfi* e 2 ^ +tt ^ * 



t 



fc l ' ~ 1 

At the surface of the body tj = 0, so that 



= -ae 



«sin§ 



2 a 2cl 



^C-j 2 +-2t L e a cos£ + e 2a ) 

(e-, 2 + e 2<t )e a cos i* 2€ 1 e 2cl 
— . g 

(C^ + cos£ * e 2a ) 



and 



V*Vtfo 



- ae^cos^ - 



> (65) 
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Therefore 



l+2e 2a U-cos -gt) 



€ X 2 + 2€ ie a cos § + ® ^+2 ^e^os £ +e 2a ) 



(64) 



The boundai-y condition at the surf aee yields Immediately 
the expression fori | ~|;] . Thus 



(65) 



where »;ns expression for / ~) 
eouations (63) „ 



given by the second of 



The expression for t ggj ^ is obtained from the equation 
for -0o Thv.s s the velocity potential can be written as s 



0 - A 1 P 1 ( l i)Q 1 (X)-> jp o (X)P o (i04-Eg(X)P 2 ( J 
Then 



Fia)Pi(ix)-t-P3a)p 3 (n 



-fllj^fA^^^+P^oj^^o)] sin £ + §P 2 (X 0 ) 



sin 2 | 



+ -^ F 3^o^ sin 3 £'- * ° {S6) 



where A x , FnO^), Fgk 0 ), ahd P 3^o) are obtained. from equations 
(59) , (60) , and (-61) » 



36 



c 

. By means of equations (63), (64), (65), and (66) 
the velocity g r of the fluid and the pressure coefficient 
P Poo 

x u3 given, respectively, by equations (50) and (49) 
can he evaluated , 
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» 

Numerical example ,- As a numerical example, consider the 

ase where *i = 0.15 and' *2 = 0,10; Then "e a - 1 +^ * ^ 
nd 



•1. 
h 




rpO 



209 ,. 5 cos _ 1540 _1 

120 8 ^ 3 317 + 75 cosZ 



> (&7) 




2& 'rp0 



5 , / , 200 \ 

s | sin *( - i «. J 



Also 



5 



1 * 



51600 



(317 + 75 cos£) 



I fcos i - 800(75.+ 317 coaJZj 
( 317 + 75 cos § ) 2 



S2 



(63) 



\2 



2a 



Tf.Q 



25 \ 
64 1 



400 



317+7-5 c 



5000 (3 3-25 cos 2£) {6g) 
(317 + 75 cos £) 2 J 



Now, 



X - 



e* + e~ a _ 41 



45 



o 2 

'i'hereforc, from equations (24) 



and ^ Q 2 -l 



9_ 
~ 40 





- 2.19723 ; 


V*o> ■ - 


1600 A u r . , 
81- ; Q o k 0 > 


5248 
' 6.5SL 




a 1.25215 ; 
= 0.82657 ; 


<V<x 0 ) » - 


18.04969 
15.99664 




Q3 <\>> 


= 0.57729 ; 


V<\>) - 


13.91684 
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and from equations (59), (60} , and (61) 

A 2 = 0.055403 ^2au); B 2 s 0.431087^ ^2au) B X = 0. 502158^2 ^2au); 
B 3 = 0.537892€ 1 2 ; F 2 (* o ) = 0.245119^ ^2au); 

Pl(\ 0 ) = 0.0791436^ (2au); Fgft 0 ) = -0.16043e€ 1 2 ^2aj) 

Then, from equation (66), 

~(^z) ~ 2aU ^0.069373 + 0.011479c ^ sin £+0.367679€ 1 sin 2$ 



2 "1 

-0.338321c j sin 3* j 



ox- with € -j — 0.15, 



~(r^l = 2au/ 0.069631 sin £ + 0^055152 sin 2£ - 0.007612 sin 3$) 

to k ~ 7 

Prom equations (65) and (68), 



(70 



V* ^/tfO L (317 f 75 cos s) 



(71) 
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Equations (75), (76), and (77) suffice to determine the pressure 
distribution over the surface of the body with the use of equations 
\*9) to (53) contained in the section entitled "Velocity and Pressure." 
Table V presents the calculated values of the various quantities 
given by equations (49) to (53), the last column containing values of the 

pressure coefficient * Figure h shows the graphs of the 

■ ipU 2 

pressure distribution for both the body of revolution (solid curve) 
and the airfoil of infinite span having the meridian profile as cross 
section (broken curve). 
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CALCULATED VALUES PROM EQUATIONS (68), (69), (70), AND (71) 



§ 

(deg) 



"55 (d|) 



1 (6 
2aU 



71=0 



0 
5 
10 
15 
20 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
130 
140 
150 
160 
170 
175 
180 



0 

-.07635 
-.15229 
-.22746 
-.30145 
-.44445 
-.57824 
-.69990 
-.80658 
-.89544 
-.96362 
-1.00813 
-1.02584 
-1.01351 
-.96807 
-.88699 
-.76911 
-.61552 
-.43043 
-.22162 
-.11164 
0 



0.30612 
.30426 
.29871 
.28953 
.27686 
.24178 
.19543 
.14049 
.08025 
.01849 
-.04072 
-.09329 
-.13552 
-.16455 
-.17889 
-.17896 
-.16735 
-.14866 
-.12888 
-.11399 
-.10988 
-.10847 



0 

-.01368 
-.02715 
-.04022 
-.05267 
-.07497 
-.09248 
-.10385 
- . 10807 
-.10469 
-.09403 
-.07724 
-.05630 
-.03379 
-.01254 
.00478 
.01615 
.02056 
.01823 
.01058 
.00548 

0 



-0.30612 
-.30426 
-.29871 
-.28953 
-.27686 
-.24178 
-.19543 
- . 14049 
-.08025 
-.01849 
.04072 
.09329 
.13552 
.16455 
.17889 
.17896 
.16735 
. 14866 
.12888 
.11399 
. 10988 
.10847 



0.09371 
.09840 
.11242 
. 13557 
. 16753 
.25600 
.37255 
.50959 
.65702 
.80216 
.93022 
1.02502 
1.07071 
1.05428 
.96917 
.81876 
.61953 
.40096 
.20188 
.06211 
.02454 
.01177 



TABLE II 

RELATIVE VELOCITY AND PRESSURE COEFFICIENTS CORRESPONDING 

TO PROFILE COORDINATES 



(deg) 
~5~ 

5 
10 
15 
20 
30 
40 
50 
60 
70 
80 
90 
100 
1-10 
120 
130 
140 
150 
160 
170 
175 
180 



x 

2a 

T.~0572 
1.0538 
1.0438 
1.0273 
1.0042 
.9390 
.8496 
.7378 
.6061 
.4573 
.2948 
.1223 
-.0556 
-.2340 
-.4074 
-.5699 
-.7149 
-.8363 
-.9279 
-.9851 
-.9997 
-1.0045 



& 



.0267 
.0530 
.0787 
.1034 
.1489 
.1872 
.2166 
.2359 
.2445 
.2425 
.2307 
.2105 
.1842 
.1540 
.1225 
.0922 
.0645 
.0403 
.0193 
.0095 



Si 

U 



.28700 
.53519 
.72700 
.86521 
1.02660 
1.09887 
1.12592 
1.12840 
1.11667 
1.09660 
1.07204 
1.04580 
1.02000 
.99609 
.97500 
.95662 
.93905 
.91741 
.84684 
.67771 

0 



1 

.91765 
.71357 
.47148 
.25141 
-.05390 
-.20751 
-.26769 
-* 27330 
-.24696 
-.20253 
-.14927 
-.09369 
-.04036 
.00781 
.04933 
.08488 
.11718 
.15836 
.28287 
.54071 

1 



TABLE III *H 
ABSCISSAS AND ORDINATES 
OF MERIDIAN PROFILE 



(percent chord) 


u 

(percent chord) 


0 . 


0 


.50 ' 


2.08 


•75 


2.53 


1.25 


3.25 


2.50 


4.62 


5.00 


6.59 


7.50 


8.06 


10.00 


9.23 - 1 


15.00 


11.04 


20.00 


12.32 


25 .00 


13.23 


30.00 


13.84 


35.00 


14.20 


40.00 


14.29 


1*5 .00 


14.12 


50.00 


13.62 


55.00 


12.69 


60.00 


11.34 


65.OO 


9.78 


70.00 


8.03 


75.00 


6.24 


80.00 


4.47 


* 85.00 


2.88 


90.00 


I.56 


95.00 


.54 


100.00 


0 
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